Deep ancestry of Bornean hunter-gatherers supports long-term local ancestry dynamics.

Borneo was a crossroad of ancient dispersals, with some of the earliest Southeast Asian human remains and rock art. The island is home to traditionally hunter-gatherer Punan communities, whose origins, whether of subsistence reversion or long-term foraging, are unclear. The connection between its past and present-day agriculturalist inhabitants, who currently speak Austronesian languages and have composite and complex genetic ancestry, is equally opaque. Here, we analyze the genetic ancestry of the northeastern Bornean Punan Batu (who still practice some mobile hunting and gathering), Tubu, and Aput. We find deep ancestry connections, with a shared Asian signal outgrouping modern and ancient Austronesian-ancestry proxies, suggesting a time depth of more than 7,500 years. They also largely lack the mainland Southeast Asian signals of agricultural Borneans, who are themselves genetically heterogeneous. Our results support long-term inhabitation of Borneo by some Punan ancestors and reveal unexpected complexity in the origins and dispersal of Austronesian-expansion-related ancestry.

The role of genetic polymorphisms of interleukin-1 (IL-1R1 and IL-1RN) in primary knee osteoarthritis in Indonesia.

This study aimed to evaluate the association of SNPs of the IL-1 family with the clinical severity of knee OA. This case‒control study was performed among 100 healthy knees and 130 osteoarthritis (OA) knees of people aged ≥ 50 years with a BMI ≥ 25 kg/m2. The possible correlations among clinical findings, radiographic evaluations, serum levels of IL-1R1 and IL-1Ra, and genotype analyses were evaluated. Three SNPs of IL-1R1, rs871659, rs3771202, and rs3917238, were associated with primary knee OA. Females with IL-1R1 SNP rs871659 allele A had a higher prevalence of primary knee OA. No correlation was found between SNPs of IL-1R1 and IL-1RN and clinical or radiologic severity or serum concentrations of IL-1R1 and IL-1Ra (p > 0.05). BMI and IL-1R1 rs3917238 genotype C/C were correlated with moderate-severe VAS scores. A correlation was also found between the EQ-5D-3L self-care dimension and obesity and between the EQ-5D-3L pain and usual activity dimensions and age ≥ 60 and obesity (p < 0.05). Radiologic severity was only associated with age ≥ 60 years (p < 0.05). We found the IL-1R1 SNPs rs871659, rs3771202, and rs3917238 to be predisposing factors for primary knee osteoarthritis. The clinical findings, radiographic severity, and serum concentrations of IL-1R1 and IL-1Ra were not correlated with these gene polymorphisms.

Human Genetic Research in Wallacea and Sahul: Recent Findings and Future Prospects.

Genomic sequence data from worldwide human populations have provided a range of novel insights into our shared ancestry and the historical migrations that have shaped our global genetic diversity. However, a comprehensive understanding of these fundamental questions has been impeded by the lack of inclusion of many Indigenous populations in genomic surveys, including those from the Wallacean archipelago (which comprises islands of present-day Indonesia located east and west of Wallace's and Lydekker's Lines, respectively) and the former continent of Sahul (which once combined New Guinea and Australia during lower sea levels in the Pleistocene). Notably, these regions have been important areas of human evolution throughout the Late Pleistocene, as documented by diverse fossil and archaeological records which attest to the regional presence of multiple hominin species prior to the arrival of anatomically modern human (AMH) migrants. In this review, we collate and discuss key findings from the past decade of population genetic and phylogeographic literature focussed on the hominin history in Wallacea and Sahul. Specifically, we examine the evidence for the timing and direction of the ancient AMH migratory movements and subsequent hominin mixing events, emphasising several novel but consistent results that have important implications for addressing these questions. Finally, we suggest potentially lucrative directions for future genetic research in this key region of human evolution.

Novel Interactions of Myristic Acid and FADS3 Variants Predict Atopic Dermatitis among Indonesian Infants.

Fatty acids exert a range of different biological activities that could be relevant in the development of atopic dermatitis (AD). This study investigated the association of glycerophospholipid fatty acids (GPL-FA) with AD, and their interactions with single nucleotide polymorphisms (SNP) of the FADS1-3 gene cluster. Among 390 infants of the Indonesian ISADI study, GPL-FA were measured in umbilical plasma (P-0y) and in buccal cells at birth (B-0y), and again in buccal cells at AD onset or one year (B-1y). Prospective and cross-sectional associations with AD were assessed by logistic regression. Interactions of GPL-FA with 14 SNP were tested assuming an additive model. AD was diagnosed in 15.4% of participants. In B-1y, C18:2n-6 was inversely associated with AD; and positive associations were observed for C18:1n-9, C20:4n-6, C22:6n-3 and C20:4n-6/C18:2n-6. There were no prospective associations with AD, however, a significant interaction between the SNP rs174449 and B-0y C14:0 (myristic acid) was observed. This study indicates that Indonesian infants with AD have increased rates of endogenous long-chain polyunsaturated fatty acid production, as well as higher C18:1n-9 levels. GPL-FA measured at birth do not predict later AD incidence; however, genotype interactions reveal novel effects of myristic acid, which are modified by a FADS3 variant.

Hypnozoite depletion in successive Plasmodium vivax relapses.

Genotyping Plasmodium vivax relapses can provide insights into hypnozoite biology. We performed targeted amplicon sequencing of 127 relapses occurring in Indonesian soldiers returning to malaria-free Java after yearlong deployment in malarious Eastern Indonesia. Hepatic carriage of multiple hypnozoite clones was evident in three-quarters of soldiers with two successive relapses, yet the majority of relapse episodes only displayed one clonal population. The number of clones detected in relapse episodes decreased over time and through successive relapses, especially in individuals who received hypnozoiticidal therapy. Interrogating the multiplicity of infection in this P. vivax relapse cohort reveals evidence of independent activation and slow depletion of hypnozoites over many months by multiple possible mechanisms, including parasite senescence and host immunity.

Episodes of Diversification and Isolation in Island Southeast Asian and Near Oceanian Male Lineages.

Island Southeast Asia (ISEA) and Oceania host one of the world's richest assemblages of human phenotypic, linguistic, and cultural diversity. Despite this, the region's male genetic lineages are globally among the last to remain unresolved. We compiled ∼9.7 Mb of Y chromosome (chrY) sequence from a diverse sample of over 380 men from this region, including 152 first reported here. The granularity of this data set allows us to fully resolve and date the regional chrY phylogeny. This new high-resolution tree confirms two main population bursts: multiple rapid diversifications following the region's initial settlement ∼50 kya, and extensive expansions <6 kya. Notably, ∼40-25 kya the deep rooting local lineages of C-M130, M-P256, and S-B254 show almost no further branching events in ISEA, New Guinea, and Australia, matching a similar pause in diversification seen in maternal mitochondrial DNA lineages. The main local lineages start diversifying ∼25 kya, at the time of the last glacial maximum. This improved chrY topology highlights localized events with important historical implications, including pre-Holocene contact between Mainland and ISEA, potential interactions between Australia and the Papuan world, and a sustained period of diversification following the flooding of the ancient Sunda and Sahul continents as the insular landscape observed today formed. The high-resolution phylogeny of the chrY presented here thus enables a detailed exploration of past isolation, interaction, and change in one of the world's least understood regions.

Interleukin-1 genetic polymorphisms in knee osteoarthritis: What do we know? A meta-analysis and systematic review.

PURPOSE: Interleukin-1 is the main proinflammatory cytokine in osteoarthritis (OA). Several single-nucleotide polymorphisms (SNPs) within the IL-1 gene cluster (IL-1ß, IL-1R1, and IL-1RN) have been determined, but their associations with knee OA remain poorly understood. The present study aimed to identify the associations between IL-1 SNPs and knee osteoarthritis. METHODS: This meta-analysis and systematic review included all comparative studies published in the MEDLINE/PubMed, Embase, Google Scholar, and Cochrane Library databases. We performed a systematic search to identify relevant studies on the evaluation of the correlation between the IL-1 gene and knee OA published up to February 2020 that met the eligibility criteria. Nine studies on a total of 2256 knees with OA and 3527 healthy knees met the eligibility criteria. Results associated with IL-1A, IL-1B, IL-1R1, and IL-1RN SNPs were extracted and compared between knees with OA and healthy knees. Methodological quality was assessed using the Newcastle-Ottawa scale (NOS). All studies with fair or good quality were included. RESULTS: The meta-analysis showed that the risk of knee OA is decreased by the IL-1RN*1 and IL-1RN*1/*1 genotypes and increased by the IL-1RN*2 and I-L1RN*1/*2 genotypes. The systematic review revealed only two studies associating the IL-1RN allele, none associating the IL-1B polymorphism, and only one study associating IL-1A and IL-1R1 polymorphisms with knee OA. CONCLUSIONS: Several IL-1RN alleles and genotypes play a role in knee OA but other genetic variations in the IL-1 region were still conflicting in its association with knee OA.

Deep ancestry of collapsing networks of nomadic hunter-gatherers in Borneo.

Theories of early cooperation in human society often draw from a small sample of ethnographic studies of surviving populations of hunter-gatherers, most of which are now sedentary. Borneo hunter-gatherers (Punan, Penan) have seldom figured in comparative research because of a decades-old controversy about whether they are the descendants of farmers who adopted a hunting and gathering way of life. In 2018 we began an ethnographic study of a group of still-nomadic hunter-gatherers who call themselves Punan Batu (Cave Punan). Our genetic analysis clearly indicates that they are very unlikely to be the descendants of neighbouring agriculturalists. They also preserve a song language that is unrelated to other languages of Borneo. Dispersed travelling groups of Punan Batu with fluid membership use message sticks to stay in contact, co-operate and share resources as they journey between rock shelters and forest camps. Message sticks were once widespread among nomadic Punan in Borneo, but have largely disappeared in sedentary Punan villages. Thus the small community of Punan Batu offers a rare glimpse of a hunting and gathering way of life that was once widespread in the forests of Borneo, where prosocial behaviour extended beyond the face-to-face community, facilitating successful collective adaptation to the diverse resources of Borneo's forests.

Mortality among healthcare workers in Indonesia during 18 months of COVID-19.

The impact of SARS-CoV-2 infections upon Indonesian health care workers (HCWs) is unknown due to the lack of systematic collection and analysis of mortality data specific to HCWs in this setting. This report details the results of a systematic compilation, abstraction and analysis of HCW fatalities in Indonesia during the first 18 months of COVID-19. HCW who passed away between March 2020 and July 2021 were identified using Pusara Digital, a community-based digital cemetery database dedicated to HCW. We calculated the mortality rates and death risk ratio of HCWs versus the general population. The analysis indicates that at least 1,545 HCWs died during the study period. Death rates among males and females HCWs were nearly equivalent (51% vs. 49%). The majority were physicians and specialists (535, 35%), nurses (428, 28%), and midwives (359, 23%). Most deaths occurred between the ages of 40 to 59 years old, with the median age being 50 years (IQR: 39-59). At least 322 deaths (21%) occurred with pre-existing conditions, including 45 pregnant women. During the first 18 months of COVID-19 in Indonesia, we estimated a minimum HCW mortality rate of 1.707 deaths per 1,000 HCWs. The provincial rates of HCW mortality ranged from 0.136 (West Sulawesi) to 5.32 HCW deaths per 1,000 HCWs (East Java). The HCW mortality rate was significantly higher than that of the general population (RR = 4.92, 95% CI 4.67-5.17). The COVID-19 pandemic in Indonesia resulted in the loss of many hundreds of HCWs, the majority of whom were senior healthcare workers. The HCW mortality rate is five times that of the general population. A national systematic surveillance of occupational mortality is urgently needed in this setting.

Genetic architecture of gene regulation in Indonesian populations identifies QTLs associated with global and local ancestries.

Lack of diversity in human genomics limits our understanding of the genetic underpinnings of complex traits, hinders precision medicine, and contributes to health disparities. To map genetic effects on gene regulation in the underrepresented Indonesian population, we have integrated genotype, gene expression, and CpG methylation data from 115 participants across three island populations that capture the major sources of genomic diversity in the region. In a comparison with European datasets, we identify eQTLs shared between Indonesia and Europe as well as population-specific eQTLs that exhibit differences in allele frequencies and/or overall expression levels between populations. By combining local ancestry and archaic introgression inference with eQTLs and methylQTLs, we identify regulatory loci driven by modern Papuan ancestry as well as introgressed Denisovan and Neanderthal variation. GWAS colocalization connects QTLs detected here to hematological traits, and further comparison with European datasets reflects the poor overall transferability of GWAS statistics across diverse populations. Our findings illustrate how population-specific genetic architecture, local ancestry, and archaic introgression drive variation in gene regulation across genetically distinct and in admixed populations and highlight the need for performing association studies on non-European populations.

Distinctive Microbiome Type Distribution in a Young Adult Balinese Cohort May Reflect Environmental Changes Associated with Modernization.

An important public health question is understanding how changes in human environments can drive changes in the gut microbiota that influence risks associated with human health and wellbeing. It is well-documented that the modernization of societies is strongly correlated with intergenerational change in the frequency of nutrition-related chronic diseases in which microbial dysbiosis is implicated. The population of Bali, Indonesia, is well-positioned to study the interconnection between a changing food environment and microbiome patterns in its early stages, because of a recent history of modernization. Here, we characterize the fecal microbiota and diet history of the young adult women in Bali, Indonesia (n = 41) in order to compare microbial patterns in this generation with those of other populations with different histories of a modern food environment (industrialized supply chain). We found strong support for two distinct fecal microbiota community types in our study cohort at similar frequency: a Prevotella-rich (Type-P) and a Bacteroides-rich (Type-B) community (p < 0.001, analysis of similarity, Wilcoxon test). Although Type-P individuals had lower alpha diversity (p < 0.001, Shannon) and higher incidence of obesity, multivariate analyses with diet data showed that community types significantly influenced associations with BMI. In a multi-country dataset (n = 257), we confirmed that microbial beta diversity across subsistent and industrial populations was significantly associated with Prevotella and Bacteroides abundance (p < 0.001, generalized additive model) and that the prevalence of community types differs between societies. The young adult Balinese microbiota was distinctive in having an equal prevalence of two community types. Collectively, our study showed that the incorporation of community types as an explanatory factor into study design or modeling improved the ability to identify microbiome associations with diet and health metrics.

Papua New Guinean Genomes Reveal the Complex Settlement of North Sahul.

The settlement of Sahul, the lost continent of Oceania, remains one of the most ancient and debated human migrations. Modern New Guineans inherited a unique genetic diversity tracing back 50,000 years, and yet there is currently no model reconstructing their past population dynamics. We generated 58 new whole-genome sequences from Papua New Guinea, filling geographical gaps in previous sampling, specifically to address alternative scenarios of the initial migration to Sahul and the settlement of New Guinea. Here, we present the first genomic models for the settlement of northeast Sahul considering one or two migrations from Wallacea. Both models fit our data set, reinforcing the idea that ancestral groups to New Guinean and Indigenous Australians split early, potentially during their migration in Wallacea where the northern route could have been favored. The earliest period of human presence in Sahul was an era of interactions and gene flow between related but already differentiated groups, from whom all modern New Guineans, Bismarck islanders, and Indigenous Australians descend. The settlement of New Guinea was probably initiated from its southeast region, where the oldest archaeological sites have been found. This was followed by two migrations into the south and north lowlands that ultimately reached the west and east highlands. We also identify ancient gene flows between populations in New Guinea, Australia, East Indonesia, and the Bismarck Archipelago, emphasizing the fact that the anthropological landscape during the early period of Sahul settlement was highly dynamic rather than the traditional view of extensive isolation.

Mitogenomes Reveal Two Major Influxes of Papuan Ancestry across Wallacea Following the Last Glacial Maximum and Austronesian Contact.

The tropical archipelago of Wallacea contains thousands of individual islands interspersed between mainland Asia and Near Oceania, and marks the location of a series of ancient oceanic voyages leading to the peopling of Sahul-i.e., the former continent that joined Australia and New Guinea at a time of lowered sea level-by 50,000 years ago. Despite the apparent deep antiquity of human presence in Wallacea, prior population history research in this region has been hampered by patchy archaeological and genetic records and is largely concentrated upon more recent history that follows the arrival of Austronesian seafarers ~3000-4000 years ago (3-4 ka). To shed light on the deeper history of Wallacea and its connections with New Guinea and Australia, we performed phylogeographic analyses on 656 whole mitogenomes from these three regions, including 186 new samples from eight Wallacean islands and three West Papuan populations. Our results point to a surprisingly dynamic population history in Wallacea, marked by two periods of extensive demographic change concentrated around the Last Glacial Maximum ~15 ka and post-Austronesian contact ~3 ka. These changes appear to have greatly diminished genetic signals informative about the original peopling of Sahul, and have important implications for our current understanding of the population history of the region.

Characterization of 29 polymorphic microsatellite markers developed by genomic screening of Sumatran rhinoceros (Dicerorhinus sumatrensis).

OBJECTIVE: The Sumatran rhinoceros is critically endangered, with fewer than 100 individuals surviving across its current range. Accurate census estimates of the remaining populations are essential for development and implementation of conservation plans. In order to enable molecular censusing, we here develop microsatellite markers with amplicon sizes of short length, appropriate for non-invasive fecal sampling. RESULTS: Due to limited sample quantity and potential lack of genome-wide diversity, Illumina sequence reads were generated from two Sumatran rhinoceros samples. Genomic screening identified reads with short tandem repeats and loci that were polymorphic within the dataset. Twenty-nine novel polymorphic microsatellite markers were characterized (A = 2.4; HO = 0.30). These were sufficient to distinguish among individuals (PID < 0.0001), and to distinguish among siblings (PID(sib) < 0.0001). Among rhinos in Indonesia, almost all markers were established as polymorphic and effective for genotyping DNA from fecal samples. Notably, the markers amplified and displayed microsatellite polymorphisms using DNA extracted from 11 fecal samples collected non-invasively from wild Sumatran rhinoceros. These microsatellite markers provide an important resource for a census and genetic studies of wild Sumatran rhinos.

Widespread Denisovan ancestry in Island Southeast Asia but no evidence of substantial super-archaic hominin admixture.

The hominin fossil record of Island Southeast Asia (ISEA) indicates that at least two endemic 'super-archaic' species-Homo luzonensis and H. floresiensis-were present around the time anatomically modern humans arrived in the region >50,000 years ago. Intriguingly, contemporary human populations across ISEA carry distinct genomic traces of ancient interbreeding events with Denisovans-a separate hominin lineage that currently lacks a fossil record in ISEA. To query this apparent disparity between fossil and genetic evidence, we performed a comprehensive search for super-archaic introgression in >400 modern human genomes, including >200 from ISEA. Our results corroborate widespread Denisovan ancestry in ISEA populations, but fail to detect any substantial super-archaic admixture signals compatible with the endemic fossil record of ISEA. We discuss the implications of our findings for the understanding of hominin history in ISEA, including future research directions that might help to unlock more details about the prehistory of the enigmatic Denisovans.

Papuan mitochondrial genomes and the settlement of Sahul.

New Guineans represent one of the oldest locally continuous populations outside Africa, harboring among the greatest linguistic and genetic diversity on the planet. Archeological and genetic evidence suggest that their ancestors reached Sahul (present day New Guinea and Australia) by at least 55,000 years ago (kya). However, little is known about this early settlement phase or subsequent dispersal and population structuring over the subsequent period of time. Here we report 379 complete Papuan mitochondrial genomes from across Papua New Guinea, which allow us to reconstruct the phylogenetic and phylogeographic history of northern Sahul. Our results support the arrival of two groups of settlers in Sahul within the same broad time window (50-65 kya), each carrying a different set of maternal lineages and settling Northern and Southern Sahul separately. Strong geographic structure in northern Sahul remains visible today, indicating limited dispersal over time despite major climatic, cultural, and historical changes. However, following a period of isolation lasting nearly 20 ky after initial settlement, environmental changes postdating the Last Glacial Maximum stimulated diversification of mtDNA lineages and greater interactions within and beyond Northern Sahul, to Southern Sahul, Wallacea and beyond. Later, in the Holocene, populations from New Guinea, in contrast to those of Australia, participated in early interactions with incoming Asian populations from Island Southeast Asia and continuing into Oceania.

Genome-wide DNA methylation and gene expression patterns reflect genetic ancestry and environmental differences across the Indonesian archipelago.

Indonesia is the world's fourth most populous country, host to striking levels of human diversity, regional patterns of admixture, and varying degrees of introgression from both Neanderthals and Denisovans. However, it has been largely excluded from the human genomics sequencing boom of the last decade. To serve as a benchmark dataset of molecular phenotypes across the region, we generated genome-wide CpG methylation and gene expression measurements in over 100 individuals from three locations that capture the major genomic and geographical axes of diversity across the Indonesian archipelago. Investigating between- and within-island differences, we find up to 10.55% of tested genes are differentially expressed between the islands of Sumba and New Guinea. Variation in gene expression is closely associated with DNA methylation, with expression levels of 9.80% of genes correlating with nearby promoter CpG methylation, and many of these genes being differentially expressed between islands. Genes identified in our differential expression and methylation analyses are enriched in pathways involved in immunity, highlighting Indonesia's tropical role as a source of infectious disease diversity and the strong selective pressures these diseases have exerted on humans. Finally, we identify robust within-island variation in DNA methylation and gene expression, likely driven by fine-scale environmental differences across sampling sites. Together, these results strongly suggest complex relationships between DNA methylation, transcription, archaic hominin introgression and immunity, all jointly shaped by the environment. This has implications for the application of genomic medicine, both in critically understudied Indonesia and globally, and will allow a better understanding of the interacting roles of genomic and environmental factors shaping molecular and complex phenotypes.

Correction to: Archaic mitochondrial DNA inserts in modern day nuclear genomes.

Following the publication of this article [1], the authors reported that the captions of Figs. 3 and 4 were published in the incorrect order, whereby they mismatch with their corresponding images.

Archaic mitochondrial DNA inserts in modern day nuclear genomes.

BACKGROUND: Traces of interbreeding of Neanderthals and Denisovans with modern humans in the form of archaic DNA have been detected in the genomes of present-day human populations outside sub-Saharan Africa. Up to now, only nuclear archaic DNA has been detected in modern humans; we therefore attempted to identify archaic mitochondrial DNA (mtDNA) residing in modern human nuclear genomes as nuclear inserts of mitochondrial DNA (NUMTs). RESULTS: We analysed 221 high-coverage genomes from Oceania and Indonesia using an approach which identifies reads that map both to the nuclear and mitochondrial DNA. We then classified reads according to the source of the mtDNA, and found one NUMT of Denisovan mtDNA origin, present in 15 analysed genomes; analysis of the flanking region suggests that this insertion is more likely to have happened in a Denisovan individual and introgressed into modern humans with the Denisovan nuclear DNA, rather than in a descendant of a Denisovan female and a modern human male. CONCLUSIONS: Here we present our pipeline for detecting introgressed NUMTs in next generation sequencing data that can be used on genomes sequenced in the future. Further discovery of such archaic NUMTs in modern humans can be used to detect interbreeding between archaic and modern humans and can reveal new insights into the nature of such interbreeding events.

Sex-linked genetic diversity originates from persistent sociocultural processes at microgeographic scales.

Population genetics has been successful at identifying the relationships between human groups and their interconnected histories. However, the link between genetic demography inferred at large scales and the individual human behaviours that ultimately generate that demography is not always clear. While anthropological and historical context are routinely presented as adjuncts in population genetic studies to help describe the past, determining how underlying patterns of human sociocultural behaviour impact genetics still remains challenging. Here, we analyse patterns of genetic variation in village-scale samples from two islands in eastern Indonesia, patrilocal Sumba and a matrilocal region of Timor. Adopting a 'process modelling' approach, we iteratively explore combinations of structurally different models as a thinking tool. We find interconnected socio-genetic interactions involving sex-biased migration, lineage-focused founder effects, and on Sumba, heritable social dominance. Strikingly, founder ideology, a cultural model derived from anthropological and archaeological studies at larger regional scales, has both its origins and impact at the scale of villages. Process modelling lets us explore these complex interactions, first by circumventing the complexity of formal inference when studying large datasets with many interacting parts, and then by explicitly testing complex anthropological hypotheses about sociocultural behaviour from a more familiar population genetic standpoint.